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We have studied theoretically the flow of superfluid 3 He- A in parallel-plate 
t> ■ geometry. The equilibrium order-parameter texture is calculated numerically 

in two spatial dimensions consisting of the coordinates along the flow direc- 
t^j- [ tion and perpendicular to the plates. The calculations have been done using 

the hydrostatic theory in the Ginzburg-Landau region and assuming a large 
external magnetic field perpendicular to the plane of calculation. We have 
studied a uniform texture and a dipole-unlocked splay soliton as initial con- 
figurations. In the former case we find the Freedericksz transition and a 
helical instability with increasing flow. In the latter case we find instability 
in the soliton. This instability is closely related to the critical velocity in the 
presence of a vortex sheet. Also, the transverse NMR frequency shift at the 
soliton has been calculated. 

O ■ PACS numbers: 67.57.Fg 
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The anisotropic A phase of superfluid 3 He exhibits unique behavior un- 
der externally applied superflow. Instead of being quantized, the circulation 
of the superfluid velocity v s in 3 He-A can have any value depending on the 
spatial variation of the order-parameter texture. The decay of superflow can 
be accomplished by generating a nonuniform texture of the orbital vector 
1, leading ultimately to the formation of continuous vortices (with no sin- 
gular core). However, energy considerations favor the uniform texture with 
constant 1 for small velocities and postpone these processes until a certain 
critical velocity is exceeded. Because of the macroscopic length scale of tex- 
ture variations and the boundary condition that fixes the direction of 1 at the 
walls, vortex formation is unaffected by extrinsic influences, such as surface 
roughness and thermal and quantum fluctuations. Therefore, in contrast 
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to superfluid 4 He, critical velocity in 3 He-A is determined by an intrinsic 
instability of the flow. 

Superflow of 3 He-A has been previously studied theoretically in many 
papers, see Refs. p] and || for references. Essentially all of them consider 
only one-dimensional flow. Here we present the first truly two-dimensional 
calculation for 3 He-A in the presence of external flow. We study a flow chan- 
nel consisting of parallel plates. With increasing flow, an initially constant 
texture becomes first modified in the Freedericksz transition, where 1 starts 
to tilt from the direction normal to the platesllfl. On increasing the veloc- 
ity further, the texture becomes unstable towards a helical deformationcEl, 
which leads to formation of vortices. However, the main point in this paper 
is to study an initially inhomogeneous texture formed by a dipole-unlocked 
solitonQ. The flow through such a domain wall was first studied by Vollhardt 
and MakS using a variational ansatz. More recently the same problem was 
studied numerically in Ref. ^. We extend these one-dimensional calculations 
to include the effect of lateral walls in a flow channel. We investigate how 
the critical velocity in the presence of the soliton depends on the width D 
of the flow channel. We find that the dependence is much stronger than for 
the helical instability in the absence of the soliton. 

The dipole-locked soliton is exceptionally interesting, because it forms 
the backbone of another topological object, the vortex sheets. Our calcu- 
lation can be used as a model for the critical velocity of the sheet, and 
it explains semi-quantitatively the measured dependence on the angular 
velocity^. Peculiar dynamical properties of the vortex sheet was observed in 
experiments^. To provide additional information to explain this behavior, 
we have also calculated numerically the NMR frequency shift which is the 
relevant experimental parameter. The shifts are determined as functions of 
flow velocity and the separation between the plates. 



2. Hydrostatic theory 

Our starting point is very similar to that in Ref. |2[ Here we repeat the 
main points. The A-phase order parameter A is a 3x3 tensor defined by two 
orthogonal unit vectors rh and n in the orbital space and a unit vector d in 
the spin space and has the formEJ (with real A) 

A^j = Adfj,(rhj + ihj). (1) 
One defines 1 = rh x n and a superfluid velocity as 
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where m.3 is the mass of a 3 He atom. The equilibrium form of the order- 
parameter texture corresponds to a minimum of the hydrostatic free energy. 
In order to avoid dealing with complicated constraints between 1 and v s , we 
express the free energy in terms of the unit- vector fields m(r), n(r), and 
d(r). For simplicity, we restrict all our calculations to the Ginzburg-Landau 
region (T c — T < T c ). The energy density consists of the dipole-dipole energy 



u-- 

the magnetic-field energy 



[d ■ m) 2 + (d • n) 2 



ft = ±A h (d • H) 2 , 



(3) 
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and the gradient energy 

ik 



** = j£?{ E [y^f + ( V *^) 2 ] + (7 - 1) [(V • m) 2 
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+ (V • n) 2 + ^(m • Vdif + J2(n • Vdi) 2 ] + 2^(Vid, 



ik 
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see Ref. 12 for details. Generally, the superfluid density is anisotropic with 
components pii and p± corresponding to the flow parallel and perpendicular 
to 1, respectively. In the Ginzburg-Landau region with the weak-coupling 
value 7 = 3, they differ by a factor of two: p± = 2py. The supercurrent 
density is given by 



Js,k 
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^2 [(7 - 2)(m k V j n j - h k Vj 



+ rhjX/ 'khj — hjVkrnj + rhjV ' jh k — hj'Vjrhk 



(6) 



A natural unit for distance is the dipole length £d = (h/2ms) wpy/Ad- Sim 



ilarly, we define units for velocity Vd = J \&/ p\\ and current jd = p\\v&. The 

free energy is minimized under the constraint that 1 must be parallel to the 
surface normal s at the walls. We restrict ourself in all following calculations 
to the large- field limit (H ^> = y Ad/Ah), where d is locked to the plane 
perpendicular to the direction of H. 

A schematic view of the situation we wish to describe is presented in 
Fig. |. In our calculations we consider a rectangular area in the xy plane 
with a length L in the direction of the flow (— L/2 < x < L/2) and width D 
in the y direction (—13/2 < y < D/2). We assume translational invariance 
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Fig. 1. The parallel-plate geometry used in the calculations. The separation 
between the plates is D and the superfluid velocity v s is in the positive x 
direction. The orientations of the external magnetic field H and the soliton 
wall (grey) are also shown. 



in the z direction. The boundary condition at the plate surfaces has the 
form rh y = h y = 0. Due to the presence of a large magnetic field H || z, 
d z = 0. 

In the numerical calculation the order parameter is defined at N x M 
discrete points spaced by Ax and Ay in x and y directions, respectively 
(NAx = L, MAy = D). The discretized free energy functional of Eqs. (||)- 
(|5|) transforms into a function of the values of rh, n and d at the lattice 
points. The minimum of this function is then sought by the following 
method: first an initial guess for the order-parameter vectors is chosen. 
Then the texture is changed towards the equilibrium distribution according 
to 

Ad = e d T d x d, 
Am = e T x rh, 

An = e r x n, (7) 

where the torques acting on the order-parameter vectors are defined as 

8F ~ 

Td = ~pr X d, 

od 

r = — xm+— xn, (8) 
r) in on 

and the (small) iteration constants ed and e D are chosen so as to achieve fast 
convergence. Also, after each iteration step the vectors are adjusted to have 
unit length and to satisfy the condition rh _L n. The iteration process is 
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repeated until the minimum of the free energy is reached (signalled by the 
vanishing of the torques). 

We generate the superflow in our model by imposing a fixed difference 
&(x = L/2) — $(x = — L/2) = A<£, which will be defined more precisely 
below. For different values of A<I>, we monitor the average induced super- 
current density 

1 f°/ 2 

(js,x) = n / dyj sx (x,y), (9) 

U J-D/2 

which is independent of x in the converged solution. We study the current- 
velocity relationship j(v), where v = (h/2m^L)A^ can be interpreted as the 
velocity of the normal component corresponding to the phase difference A$ 
(Ref. H). The critical velocity v c of the flow instability is determined from 
the condition 

Y V U=V c =0. (10) 



J 



3. Uniform initial state 

Before studying the soliton case, it is necessary to understand the flow 
response in the absence of the soliton. We begin by considering the simplest 
order-parameter structure in 3 He-A between parallel plates. For the case of 
zero velocity, v = 0, the free energy is minimized by a uniform texture, where 
1 || d || y. This homogeneous configuration is a simultaneous minimum of 
/gi /di and /h, and also satisfies the requirement 1 || s at the boundaries. 
Next we introduce a flow in the system by a requirement 

A(x = L/2,y) = e iA *A(x = -L/2,y) 

= e 2ivmsL/h A (x = -L/2,y). (11) 

In the presence of flow, 1 tends to turn towards the flow direction because 
the component of the superfluid-density tensor along 1, py , is smaller than 
the perpendicular one, p±. This would require the formation of an inhomo- 
geneous texture, because at the surfaces of the plates 1 is rigidly anchored 
perpendicular to the flow. Therefore, for small enough velocities the uni- 
form texture 1 1 v s persists as the equilibrium configuration, with m and 
n turning around constant 1 to achieve the required phase difference or, 
equivalently, A= const, x exp(2ivm^x/h). 

The transition where 1 starts to deflect from the plate normal y is known 
as Freedericksz transition, and v-pt denotes the corresponding velocity. The 
Freedericksz transition has beenstudied extensively in the case where it is 
induced by magnetic field H || yEJ. In the present case, it was first calculated 
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Fig. 2. The variation of 1 (solid arrows) and d (dashed arrows) for D = 6£d> 
when vyt ~ 0.51v^. Both fields are translationally invariant in x direction. 
The textures are for a) v < vf t (homogeneous texture), b) v = 0.55vd, and 
c) v = 1.3ud- 



by deGennes and RainerB, who neglected the dipole-dipole interaction (|3|) 
and obtained 

^Fr = ^, (12) 

where c = (irh/2m3)y/3/4. This result holds in the limit D but in a 

wider slab there is additional textural rigidity because d tends to follow 1. 
This has been calculated both exactly and by variational ansatz by FetterQ. 
In the limit D 3> £d there is complete dipole-locking and one obtains (|12|) 
with c = (irh/2m3) At intermediate D's VFr monotonically interpo- 
lates between these limits (see Fig. || below). 

Typical 1 and d textures for different values of v are presented in Fig. |2|. 
On increasing the velocity further beyond vp T , 1 bends more and more until 
most of the texture (excluding the vicinity of the surfaces) has aligned itself 
with the flow direction. Finally, at v = v^, the texture becomes unstable 
towards helical disturbances that break the translational invariance in the 
x coordinated, see Fig. ||. We have not made systematic studies of this 
transition here, but it seems that there is no substantial difference to the 
corresponding one-dimensional calculationa, which corresponds to the limit 
D — > oo. The critical velocity increases slightly because of the presence of 
the lateral walls, 

v h (D) w ^h(oo) + (13) 

where Wh(oo) = 1.26t>d (Ref. [|) and a ~ 3v^^. Another question is the 
stability of these helical textures. It has been found in one dimension that 
the helical textures are unstable towards nucleating vortices in the vicinity 
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Fig. 3. Helical disturbances of 1 (solid arrows) and d (dashed arrows) for 
D = at v = 1.3ud- A length of L = 15£d in the direction of flow is 
presented, comprising three wavelengths of the helix. 



of r c &EI. Whether this is true also in our restricted geometry remains open. 

A useful way of expressing the response an externally applied flow is 
through the cur rent- velocity diagram j(v). A typical example of such a 
diagram is presented in Fig. ^. At low velocities the uniform texture with 
11 v s can be seen as a linear slope j = p±v. At the Freedericksz transition 
part of the 1 texture starts to deflect towards v s , causing a decrease in the 
slope of j. With increasing velocity the texture tends to align with v s and 
the slope approaches p\\. 



4. Flow with a soliton 

The complicated order-parameter structure of 3 He-A allows for a va- 
riety of different topological defects. One variety of these are solitons, i.e. 
two-dimensional domain-wall-like structures separating two energetically de- 
generate minima of the dipole-dipole interaction (0). On one side of the 
soliton 1 and d are parallel to each other, on the other side antiparallel. In 
our geometry the soliton has the so-called splay structure, i.e. the magnetic 
field is in the plane of the soliton, see Fig. || Because of the orienting effect 
of the walls, far away from the soliton a uniform 1 = y texture is approached, 
with d = — y on the left side and d = y on the right side of the wall. 



The boundary condition (11) we used to generate flow in the uniform 
case is not consistent with the presence of a soliton. Instead, it is possible to 
show that the combined operation of translation by L in the flow direction 
and reflection in the y coordinate is a symmetry operation for the soliton 
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v ( v d ) 

Fig. 4. The cur rent- velocity diagram for D = 8£d- The Freedericksz 
transition shows up as a kink at v ~ OAv^ and the curve ends to the onset 
of helical texture at v ~ lSv^. The lines j = p±v and j = puv are shown as 
guides to the eye. 



texture. Therefore, a flow velocity equal to v can be achieved by requiring 
A(x = L/2, y) = e 2ivmaL / H a y A (x = -L/2, -y)a y , (14) 

where a y is a diagonal matrix with elements 1,-1, and 1. Otherwise the nu- 
merical calculations follow the same procedure as in the case with a uniform 
initial texture. 

In order to understand the response of the soliton wall to an external 
flow, it is advantageous to first study a simpler case in one dimension x along 
the flow in the bulk. We consider the range from x = —L/2 to x = L/2, 
over which the direction of 1 is reversed from 1 = — x to 1 = x, see Fig. [6|. 
In the absence of a magnetic field, it is possible to show that, starting with 
a texture with zero phase difference over the length L, a texture with a 
phase difference equal to A<I> can be formed by simply rotating the initial 
texture an angle A<3?/2 around x. Furthermore, also the final texture carries 
zero supercurrent. The situation changes if a magnetic field H _L x is intro- 
duced, because it prevents 1 from rotating freely through the dipole-dipole 
interaction. As a result, a texture corresponding to a nonzero supercurrent 
is formed. When the phase difference exceeds a critical value, the texture 
undergoes an abrupt rotation through a large angle, again leading to vortex 
creation somewhere on the sides. Roughly the same scenario takes place also 
in the two-dimensional parallel-plate geometry, with some modifications. Far 
away from the soliton plane and at the solid walls the texture is locked to 
a uniform 1 = y configuration and cannot rotate when the flow is applied. 
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Fig. 5. The variation of 1 (solid arrows) and d (dashed arrows) in a splay 
soliton structure for D = 8£d and zero velocity. On the left side 1 = — d and 
on the right side 1 = d. 



This results in an additional rigidity opposing the rotation of the soliton. 

The calculation of the critical velocity for the soliton presents some 
problems because, due to the breakdown of translational invariance in the 
x direction, the flow response of the system depends on the length L of our 
computational region. We define the critical velocity as 

v c = lim (h/2m 3 ) A ®f L \ (15) 

L^OO Li 

where A<£ c corresponds to the phase difference for which the maximum cur- 
rent j c (L) is achieved. A typical form of the current-velocity diagram at 




Fig. 6. Flow response of a one-dimensional soliton structure (schematic). 
A phase difference of A<£ can be accomplished by rotating texture an angle 
of A<J>/2 around the flow direction. 



J. Kopu and E. V. Thuneberg 




v(v d ) 

Fig. 7. The cur rent- velocity diagram for a splay soliton texture for D = 
L = 8£d (solid line). The diagram in the absence of the soliton for a channel 
of the same width is also shown (dashed line). 



finite L is shown in Fig. |7| With increasing L (thus increasing the relative 
length of the soliton-free region) the current-velocity curve approaches the 
translationally invariant behavior discussed in the previous chapter. How- 
ever, it turns out that v c (L) shows appreciable dependence on L even for 
L ^> On the other hand, the value of the critical current j c (L) only varies 
on a much smaller scale (on the order of a few percent). The dependences 
are presented in Fig. |8|. Therefore, in order to avoid lengthy calculations 
involving a huge number of lattice points, we calculate the critical current 
using a moderate-sized lattice and read the corresponding velocity v c from 
the current- velocity curve without the soliton (such as in Fig. ||). 

The critical current j c for the splay soliton texture is presented in Figs. ^ 
and 10 as a function of the plate separation D. In Fig. || we have plotted 
the product j c D. We have also included for comparison the numerically 
calculated jFr-D and the analytic limits (^) for the Freedericksz transition. 
In order to have a more accurate look at the limit of large plate separations, 
we present j c and 2v c as functions of inverse D in Fig. [u]. As long as the 
critical currents exceed jp r = 2vp T (i.e. when D < 12£d)i the dimensionless 
critical currents and velocities are simply related through a factor of two: 
Jc/jd = 2w c /wd- However, for large plate separations j exceeds jF r , and the 
critical velocity decreases more slowly (see Fig. Q). 

From the figures we find that, for intermediate values of D, j c obeys 
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Fig. 8. The critical velocity u c and current j c for a splay soliton as functions 
of L for D = 4£ d . 



reasonably well the relation 

jc « £ + 6, (16) 

with o = 4.5j d £ d , 6 = 0.17j d for D < 12£ d , and a = 5.1j d £ d , & = 0.11j d 
for D > 12£ d . Fig. [h] suggests that t> c levels off at the limit of large D to 
0.15f d < v^ nlk < 0.2u d , in agreement with our previous result obtained from 
a one-dimensional treatment of a dipole-unlocked solitona. Note also that 
v c/vd = jc/jd in the limit D — > oo. The main conclusion from the calcula- 
tions above is that the critical velocity in the presence of a soliton is small 
and depends essentially on the width D of the channel. The dependence is 
much stronger than for the helical instability in the absence of the soliton. 

The calculation above gives a qualitative explanation for the measured 
critical velocity of a vortex sheet. The vortex sheet consists of a soliton 
backbone to which continuous vorticity is addedtl In a rotating container 
the sheet has end lines at the side walls of the container. During angular 
acceleration there is counterflow through the vortex-free end pieces of the 
sheet, and the sheet can grow only when new vorticity is nucleated there. It 
is found experimentally that the critical velocity for the growth of the vortex 
sheet depends essentially on the angular velocity VL, increasing with increas- 
ing This can now be understood because the width of the vortex-free 
soliton pieces decreases with increasing fi, which corresponds to decreasing 
D in a channel and thus increasing t; c . The channel model is still a rather 
crude approximation for the end piece of a vortex sheet, but in spite of that it 
gives the right order of magnitude for the critical current and its dependence 
on Q. 
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Fig. 9. The critical current j c times the plate separation D for a splay 
soliton as a function of D (solid line) Also shown is j-ptD corresponding to 
the Freedericksz transition of a uniform texture (dashed line) , together with 



analytic results (12) for the limiting cases of dipole-locking (upper dash- 
dotted line) and for d = constant (lower dash-dotted line). 
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Fig. 10. The critical current j c and twice the critical velocity 2v c for a splay 
soliton as a function of 1/D. 
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We can argue analogously for the continuous vortex lines. Since we 
found only weak dependence of the helical instability on the channel width 
(]j~3|), the critical velocity for creating vortex lines should be essentially inde- 
pendent on the rotation velocity. This is what is found experimentally0. 



5. Nuclear magnetic resonance 

NMR techniques are effective in obtaining information about the order- 
parameter distribution in superfluid 3 He. The frequency of the transverse 
NMR absorption peak is given by0'SJ§ 

J 1 = loI + R\u\, (17) 

where lvq is the Larmor frequency, the A-phase longitudinal resonance 
frequency, and R\ is determined by solving the eigenvalue problemD-3 

- (2V 2 - V • 11 • V)^ + Uip = (R\ - l)ip (18) 

for the transverse fluctuations of d given by the wave function if). The 
potential for the fluctuations is given by 

U = -P z -(lx d) 2 - 2(W) 2 + (1 • V0) 2 , (19) 

where d = x cos 9 + y sin 9 corresponds to the unperturbed distribution in 
the presence of a large magnetic field H || z. An approximation for the 
lowest eigenvalue and the corresponding bound state of fluctuation induced 
by the rf field can be obtained by the following variational principle: 

p2 -i . /d 3 r(2|vyf-|i-vyf + £/|yf) 

R ± -l = mm ^ . (20) 

i> J d d r\^\ 2 

For a uniform state with d || 1 _L H, the potential in Eq. ( |l9l) vanishes, tp is 
constant and the absorption occurs at the bulk frequency given by R\ = 1. 
However, dipole-unlocked regions such as the splay soliton define a potential 
well, and usually give rise to the presence of bound states with frequencies 
corresponding to R\ < 1. The effective volume of absorption associated 
with a given mode in a homogeneous rf field is given by 

| fd 3 n/>(r)expHfl(r)]| 2 
= • (21) 

In the lowest mode 2 is expected to accumulate near the minima of the 
potential energy U, which in our problem occur at the two points where the 
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Fig. 11. Transverse NMR frequency shift R\ as a function of the current j 
for two different plate separations of 6£d (solid line) and 10£d (dashed line). 



soliton joins the plates. Due to the broken reflection symmetry in y — > —y 
(s • 1 changes sign at the wall, see Fig. |5|) two effects arise: the potential wells 
have different depths, and the directional derivatives |1 • V?/>| 2 of Eq. ( p0[ ) 
give different contributions near the two joining points. As a result, one of 
the points is singled out as D is increased. 

Recent experiments on vortex sheet dynamics have exhibited unex- 
pected behavioral. A possible explanation involves a rotating state that 
consists of several sheets and multiple contact points with the container. 
The sheets are joined to the walls by vorticity-free soliton pieces. An im- 
portant question arising in the explanation is whether these pieces have the 
same resonance frequencies than the sheet itself. Figure 11 shows R 2 ^ as 
a function of the flow velocity for two different plate separations. The fre- 
quencies we obtained are considerably lower than those associated with bulk 
splay soliton, calculated e.g. in Ref. [| for zero velocity as R\ = 0.677. In 
fact, the frequencies are rather close to the ones expected for a vortex sheet, 
previously calculated to be Rj_ = 0.46-0.48, depending on the models! 

In addition to the frequency shift associated to a given mode, another 
important quantity is the relative intensity of the mode. This can be ex- 
pressed in terms of the effective volume of absorption, Vnmr, as defined 
in Eq. (]2l]). For small values of D, \ip\ 2 is distributed almost evenly along 
the entire soliton, and Vnmr grows linearly with increasing D. After D is 
increased enough, \ip\ 2 concentrates at the joining point of the soliton and 
one of the plates with y = — D/2 (s • 1 = 1), resulting in a constant Vnmr- 
Figure ^ shows the effective area of absorption (volume per unit length in 
z direction) for zero phase difference, v = 0. For comparison, we have also 
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Fig. 12. The effective area of absorption Snmr °f the lowest eigenmode as a 
function of the plate separation D for zero velocity (solid line). Also shown 
is the area of absorption for a bulk splay soliton of length D (dashed line). 



drawn the area of absorption for a bulk si 
to have a absorption thickness of 7.3£a" 



ay soliton of length D, calculated 
. The absorption thickness of a 



vortex sheet is approximately 6-8£d depending on the angular velocity^ 



6. Conclusions 

We have studied the flow stability of superfluid 3 He-A in the vicinity 
of solid walls by calculating numerically the order-parameter distribution in 
two spatial dimensions. The emphasis of this work was to study the flow in 
the presence of a splay composite soliton. We have determined the critical 
current and velocity for the creation of continuous vortices as functions of 
the separation between the walls. The results give a qualitative explanation 
to the experimentally measurable critical velocity for the growth of a vortex 
sheet. The transverse NMR frequency shift B?\_ for the lowest eigenmode of 
the dipole-unlocked soliton was found to be rather close to that of a vortex 
sheet. 
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